home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / sptsv.z / sptsv
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. SSSSPPPPTTTTSSSSVVVV((((3333FFFF))))                                                            SSSSPPPPTTTTSSSSVVVV((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SPTSV - compute the solution to a real system of linear equations A*X =
  10.      B, where A is an N-by-N symmetric positive definite tridiagonal matrix,
  11.      and X and B are N-by-NRHS matrices
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE SPTSV( N, NRHS, D, E, B, LDB, INFO )
  15.  
  16.          INTEGER       INFO, LDB, N, NRHS
  17.  
  18.          REAL          B( LDB, * ), D( * ), E( * )
  19.  
  20. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  21.      SPTSV computes the solution to a real system of linear equations A*X = B,
  22.      where A is an N-by-N symmetric positive definite tridiagonal matrix, and
  23.      X and B are N-by-NRHS matrices.
  24.  
  25.      A is factored as A = L*D*L**T, and the factored form of A is then used to
  26.      solve the system of equations.
  27.  
  28.  
  29. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  30.      N       (input) INTEGER
  31.              The order of the matrix A.  N >= 0.
  32.  
  33.      NRHS    (input) INTEGER
  34.              The number of right hand sides, i.e., the number of columns of
  35.              the matrix B.  NRHS >= 0.
  36.  
  37.      D       (input/output) REAL array, dimension (N)
  38.              On entry, the n diagonal elements of the tridiagonal matrix A.
  39.              On exit, the n diagonal elements of the diagonal matrix D from
  40.              the factorization A = L*D*L**T.
  41.  
  42.      E       (input/output) REAL array, dimension (N-1)
  43.              On entry, the (n-1) subdiagonal elements of the tridiagonal
  44.              matrix A.  On exit, the (n-1) subdiagonal elements of the unit
  45.              bidiagonal factor L from the L*D*L**T factorization of A.  (E can
  46.              also be regarded as the superdiagonal of the unit bidiagonal
  47.              factor U from the U**T*D*U factorization of A.)
  48.  
  49.      B       (input/output) REAL array, dimension (LDB,N)
  50.              On entry, the N-by-NRHS right hand side matrix B.  On exit, if
  51.              INFO = 0, the N-by-NRHS solution matrix X.
  52.  
  53.      LDB     (input) INTEGER
  54.              The leading dimension of the array B.  LDB >= max(1,N).
  55.  
  56.      INFO    (output) INTEGER
  57.              = 0:  successful exit
  58.              < 0:  if INFO = -i, the i-th argument had an illegal value
  59.              > 0:  if INFO = i, the leading minor of order i is not positive
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. SSSSPPPPTTTTSSSSVVVV((((3333FFFF))))                                                            SSSSPPPPTTTTSSSSVVVV((((3333FFFF))))
  71.  
  72.  
  73.  
  74.              definite, and the solution has not been computed.  The
  75.              factorization has not been completed unless i = N.
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.